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摘要 : ”首先 根据 传统 节点 形 函 数 定义 域 的 范围 ， 提 出 分 节点 和 全 节点 的 概念 ， 基于 传统 三 
角形 板 壳 单元 分 节点 形 函 数 , 通过 采用 平行 四 边 形 方 法 , 对 分 节点 形 函数 的 定义 域 关于 坐标 
原点 进行 的 延 拓 ， 构 建 了 全 节点 形 函 数 ， 然后 ， 将 此 全 节点 形 函 数 伸缩 后 在 单元 区 域 上 进 
行 平移 ， 生 成 基 函 数 系 ， 即 多 分 辨 分 析 的 简单 而 明晰 的 数学 基 ， 同 时 ， 该 基 张 成 了 位 移 子 
空间 序列 ， 由 此 建立 了 带 有 分 辨 率 常 数 的 理性 多 辨 分 析 的 概念 , 并 在 此 基础 上 ， 进 一 步 构建 
了 多 分 辨 率 三 角形 板 壳 单元 及 其 单元 法 ， 最 后 ， 通 过 算 例 验证 ， 可 以 得 到 如 下 结论 , 分辨 率 
定义 了 区 域内 节点 的 等 距 均匀 分 布 形 式 ; 传统 三 角形 板 壳 单元 为 单 分 辨 率 单 元 ， 为 本 文 单元 
的 一 个 特例 ; 由 于 域内 为 全 节点 ， 所 以 本 文 建立 的 结构 计算 模型 是 整体 化 模型 而 不 是 传统 
的 网 格 碎 片 化 模型 ; 本 文 方法 为 理性 多 辨 分 析 方 法 ， 其 计算 效率 高 于 其 他 非 理性 多 辨 分 析 方 
法 且 可 以 完全 统一 其 他 的 非 理性 多 辨 分 析 方 法 ; 结构 数值 分 析 清 晰 度 的 高 低 是 由 分 辨 率 的 大 
^v x BJ fi aF 0416 Bü 96 s 连续 的 全 节点 形 函数 为 传统 的 节点 总 刚 可 由 共同 节点 的 单元 刚度 通 
过 人 工 闭 加 组 合 方式 而 得 到 的 这 种 处 理 方法 提供 了 理论 依据 。 

X 键 dd. 三 角形 板 壳 单元 ; 分 节点 ; 全 节点 ; 位 移 子 空间 序列 ; 理性 多 辨 分 析 ; 分 
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‘= 有 限 元 方法 [1] 已 成 为 当今 各 个 领域 进行 科学 研究 和 工程 计算 普遍 采用 的 通用 方法 ， 它 


对 社会 生产 各 个 行业 产生 了 深远 的 影响 。 板 吝 单 元 是 有 限 元 法 中 一 种 常用 的 单元 类 型 , 其 用 
来 计算 工程 中 经 常 采用 的 板 壳 结构 。 常 用 的 板 壳 单元 一 般 可 分 为 矩形 和 三 角形 单元 , H-T 3E 
形 单元 不 适合 边界 几何 形状 比较 复杂 的 结构 ， 因 此 , 在 实际 应 用 中 , 三 角形 单元 得 到 了 更 为 
广泛 的 应 用 , 三 角形 单元 在 复杂 网 格 划分 时 ， 可 采用 计算 几何 学 中 的 Delaunay 三 角形 算法 ， 
其 提高 了 网 格 划分 效率 ; 另外 ,在 需要 网 格局 部 加 密 时 ， 三 角形 单元 也 更 显 便捷 ， 尽 管 存在 
这 样 的 优势 ,但 由 于 没有 一 个 可 以 遵循 的 法 则 ， 可 以 预先 确定 三 角形 网 格 的 尺寸 大 小 、 各 个 
节点 数量 和 位 置 ， 因此 ， 对 于 工程 中 的 有 具体 问题 ， 即 使 是 经 验 丰富 的 技术 人 员 ， 要 选择 合 
适 的 网 格 方案 也 很 困难 ， 往 往 需 要 多 次 反复 才能 达到 指定 的 求解 精度 。 
三 角形 网 格 节点 数量 和 位 置 的 确定 无 法 则 可 循 的 问题 来 自 于 有 限 元 理论 自身 所 固有 的 
缺陷 ， 其 主要 表现 在 两 个 方面 : 一 是 三 角形 单元 节点 形 函数 的 定义 域 限定 在 单元 区 域内 ， 导 
臻 单元 各 完整 节点 (全 节点 ， 如 图 2 中 的 涂 黑 原点 〉 被 分 割 成 分 离 节 点 分 节点 ， 如 图 1b 中 
涂 黑 部 分 ),， 所 以 , 单元 网 格 划分 的 过 程 就 是 将 整个 结构 体 从 各 单元 节点 处 进行 离散 的 过 程 ; 
其 二 是 不 确定 性 源 自 于 节点 分 布 的 任意 性 , 即 网 格 划分 可 采用 节点 规则 均匀 分 布 的 形式 , 也 
可 采用 节点 随机 不 规则 分 布 的 形式 ， 终 其 原因 是 缺少 一 种 控制 节点 数量 及 其 分 布 形式 的 法 
则 ， 即 没有 一 个 坚固 数学 基 的 支撑 。 一 般 而 言 ， 多 关 分 析 是 指 曝光 的 信息 细节 (节点 ) 数量 


| 
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可 以 某 种 方式 任意 改变 的 技术 ， 以 均匀 分 布 的 方式 (有 坚固 数学 基 支 撑 ) 而 改变 节点 数量 的 
多 辨 分 析 可 称 为 理性 多 辨 分 析 ; 以 随机 任意 不 规则 分 布 方 式 (无 数学 基 支 撑 ) 而 改变 节点 数 
量 的 称 为 非 理 性 多 辨 分 析 。 显 然 ， 三 角形 网 格 中 完整 信息 细节 (全 节点 ) 被 撕 裂 分 离 以 及 非 


理性 多 辨 分 析 是 有 限 元 理论 两 个 固有 的 漏洞 。 


在 过 去 的 三 十 余年 中 , 为 了 消除 网 格 划分 带 来 的 问题 , 各 国学 者 相继 提出 了 各 种 办 新 的 


方法 。Belytschko 教 授 等 在 1994 年 提出 了 无 网 格 法 [2][3]， 该 法 采用 目标 点 周围 
数 的 条 件 来 确定 多 项 式 系数 的 方法 , 由 此 构建 此 目标 全 节点 的 形 函 数 ; Ko 教授 等 用 
用 Daubechies 小 波 尺度 函数 矢量 积 的 方法 建立 了 域内 全 节点 形 函 数 ， 从 而 提出 了 小 波 有 限 元 
法 [4]; Hughes 教 授 等 于 2005 年 建立 了 等 几何 法 [5][6]， 该 法 采用 样 条 函数 的 矢量 积 
建 了 域内 全 节点 形 函 数 ， 可 以 看 出 ， 这 些 方法 有 一 个 共同 的 特点 ， 就 是 构建 了 域内 全 节点 形 
函数 ， 从 而 建立 了 结构 的 整体 计算 模型 ， 而 不 是 传统 有 限 元 需要 划分 网 格 的 离散 计算 模型 ， 


点 应 满足 形 函 
E1995 年 利 


的 方法 构 


所 以 , 采用 这 些 方法 时 无 需 划 分 网 格 , 因此 它们 在 解决 某 类 特殊 工程 问题 时 能 够 显示 出 它们 
各 自 的 优势 ， 但 是 ， 由 于 此 类 方法 在 节点 布置 方式 上 比较 随意 ， 没 有 数学 基文 撑 ， 所 以 ， 它 


们 的 计算 依旧 属于 非 理性 多 辨 分 析 。 


针对 节点 布置 方式 上 随意 , 没有 一 个 数学 基 支 撑 的 问题 ， 本 研究 提出 
KE， 该 法 通过 延 拓 分 节点 形 函 数 ， 形 成 基本 全 节点 形 函 数 ， 并 在 此 基础 上 ， 采 用 4 


YE 


了 多 分 辨 率 单元 方 


]4 EE 


函数 在 单元 区 域 上 进行 收缩 和 平移 的 方法 ， 构 造 了 一 个 基 函 数 系 ， 也 就 


是 一 个 简洁 | 


的 数学 


基 ， 从 而 构建 了 一 个 位 移 子 空间 序列 , 并 由 此 形成 了 一 个 细 新 的 带 有 分 辨 率 常数 〈 应 ) 的 多 


辨 分 析 ， 同 时 ,分辨 率 定 义 了 节点 均匀 分 布 的 法 则 ， 它 可 以 通过 调节 单元 分 辨 率 常 数 来 自由 


调整 节点 的 数量 ,一旦 确定 了 分 辨 紊 常数， 那么 节点 的 数量 和 位 置 也 就 相应 给 定 了 ， 从 而 解 


决 了 传统 有 限 元 和 一 些 其 它 方法 中 所 存在 的 节点 数量 及 其 位 置 不 确定 的 问题 。 分 辩 率 是 指 能 


够 识别 出 两 点 信息 细节 《物理 量 值 ) 存在 差别 时 ， 该 两 点 之 间 的 最 小 距离 。 因 此 ， 统 一 的 分 
状 率 要 求 区 域内 的 相 邻 节点 之 间 应 保持 相同 的 距离 。 由 此 可 以 看 出 , 多 分 辩 率 单元 法 是 具有 


坚固 数学 基文 撑 的 理性 多 辨 分 析 方 法 , 无 网 格 , 它 运 用 分 辨 率 常 数 和 分 辨 率 来 而 


外 定 节 点 的 数 


量 及 其 位 置 , 即 有 规 可 循 。 本 文 就 三 角形 薄板 这 单元 进行 多 分 辨 分 析 计 算 公 式 的 推导 ， 为 解 


决 工程 中 板 壳 结 构 的 计算 问题 提供 另 一 种 强 有 力 的 工具 。 


1 基本 全 节点 形 函数 


有 一 三 角形 薄板 壳 单 元 , 建立 如 图 1a 所 示 的 直角 坐标 系 , 薄板 这 单元 的 中 面 直角 


如 图 15 所 示 。 


坐标 系 


+ 


图 1. 三 角形 薄板 壳 单 元 


Fig 1. A triangular shell element 
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由 此 ， 我 们 可 以 写 出 如 图 坐标 系 下 三 角形 三 条 边 
对 于 底线 12 


2-0 
h 


对 于 边线 23 


式 中 a,b 分 别 为 水 平 轴 和 竖 轴 的 截 距 . 
对 于 边线 13 


那么 ， 三 角形 板 壳 单 元 上 任意 一 点 的 位 移 可 写 为 : 


HEN +Y N 0 
i=l 


But. vt. 


2. 3) 沿 坐 标 轴 x、y、z 方 向 的 位 移 ;，0.,、 


单元 节点 ;处 的 平面 
单元 节点 1 处 绕 x、y 轴 的 转角 形 函 数 。 
Est OD 中 的 各 个 形 函 数 可 写 为 [1]: 
NP =L, 
NP =L+CL+GL-LÉ-LE 


ns 


1 1 
N° =—b, (aL “Ë ^ LLL +b, (zz + 5 LLL 


N’ = - (£t egets ee [nent] 


b-y-y, G =X,-X, b =y; -y> =X 


HHL. L,. 、 万 分别 是 三 角形 


由 式 (1) 、 (2). (3) 可 得 到 如 下 关系 式 


Ww 分别 为 沿 坐 标 轴 x、y、z 方 向 的 位 移 ;,、 


的 无 量 纲 


薄膜 内 位 移 形 函 数 ;，N? 为 单元 节点 7 处 的 挠 | 


=X 


形 单元 三 个 节点 的 面积 坐标 。 


(1,2,3 


函数 方程 ; 


0, 分别 为 单元 节点 1 处 绕 t、 


UE PR s N? ` 
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(1) 


(2) 


(3) 


(4) 


的 转角 ，N” 为 


b ` 
N* 分 别 为 
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x y 1 1 y 
L=1-| —+=|, cus ee , L,=—, xyeD (5a. b.c) 
1 E 3 2 h 3 șa y 1 
显然 有 +L+L=1， 且 在 其 三 角形 单元 定义 域 D A, IL. L, 1,20 


接 下 来 ,我 们 可 以 采用 上 面 的 三 角形 板 壳 单元 分 节点 形 函 数 ,来 构建 全 节点 形 函 数 。 如 
图 lb 所 示 ， 三 角形 单元 的 节点 为 单元 定义 域 〈 阴 影 部 分 ) 所 分 割 〈 涂 黑 部 分 )， 即 单元 分 裂 
了 单元 节点 ， 因 此 ,传统 有 限 元 中 的 单元 网 格 划分 是 从 各 个 节点 处 离散 了 结构 整体 。 为 了 建 


立 非 离 散 的 整体 结构 计算 模型 ， 必 须 构 建 完整 节点 -全 节点 形 函 数 ， 此 函数 的 构建 关键 是 延 
拓 节 点 形 函 数 的 定义 域 ， 即 定义 域 应 闭 包 该 节点 ， 因 此 ， 针 对 三 角形 单元 , 我 们 可 以 采用 以 
原 三 角形 为 基础 ， 围 绕 原 点 (节点 ) 作 一 系列 的 平行 四 边 形 ， 即 采用 登 加 相同 形状 和 大 小 的 


三 角形 方式 扩展 此 节点 的 定义 域 ， 最 终 得 到 如 图 2 中 的 六 边 形 ， 显 然 ， 图 中 原点 处 的 节点 为 
其 定义 域 〈《 阴 影 部 分 》 闭 包 ， 此 节点 应 为 全 节点 ， 下 面 给 出 该 基本 全 节点 的 形 函 数 。 
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图 2. 三 角形 定义 域 的 延 拓 


Fig 2. Extension of a triangular domain 
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RPN”, N’, N5. NO 分 别 为 传统 三 角形 板 壳 单元 分 节点 i (二 1、2、3》 处 的 形 函 数 ， 
即 基本 全 节点 域 由 六 个 三 角形 分 节点 域 所 县 加 而 成 。 其 中 有 
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到 此 ， 得 到 了 基本 全 节点 形 函 数 ， 其 形状 如 图 3 所 示 


201801.00046v4 


chinaXiv 


a. @(x, y) b. g(x,y) 


图 3. 基本 全 节点 形 函 数 
Fig 3. Plot of the basic shape function 


显然 ， 该 基本 全 节点 形 函 数 是 连续 的 且 具 有 Kronecker delta 特性 : 
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2 位 移 子 空间 序列 


为 了 实现 多 分 辩 率 的 计算 分 析 , 我 们 须 构建 


壳 单 元 ， 当 单元 分 辨 率 常 数 RL(resolution level) s(n + 1)(n + 2) Bf, W^ 可 取 为 


€ [I9 us Oe | (11) 


n,rs n,nn 


单元 的 多 分 辨 率 位 移 子 空间 序列 。 通 过 函数 
的 伸缩 、 平 移 可 以 组 建 子 空 间 序 列 的 基底 函数 系 ,由 此 ,将 基本 全 节点 形 函 数 在 区 域内 进行 


伸缩 、 平 移 ， 由 此 形成 单元 的 多 分 辨 率 位 移 子 空间 序列 的 基 函 数 向 量 亚 , ， 对 于 三 角形 薄板 
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Q, ,7 P(X — rast an-s) š doe (r7 52 Jam = sh š 


ds meros mro. ss 76 0 (r7 52 Jom. n 为 正 整数 ， 表 


示 三 角形 各 个 边 的 伸缩 倍数 ; r= 0,1,2,3…n,s = 0,1,2,3.…n 表示 节点 位 置 平移 的 参数 ; 


这 里 n>r>s, (n-ra)e[-aa] (ny-sh)e[-h,h] .x;y € D,. 


WA faeere -2,2 |x| 2.2 cae, amana 
nn 


nn 


雍和 到 元 区域 上 的 所 有 点 (o2) Sa Je 由 此 ， 即 可 构造 出 单元 区 域 D. 
n n 


内 各 个 全 节点 的 形 函数 。 

如 图 4 所 示 , 对 于 三 角形 单元 域内 的 全 节点 形 函 数 , 可 由 基本 全 点 形 函数 经 过 收缩 和 移动 
来 得 到 ; 对 于 三 角形 单元 域 边界 节点 (B) 的 形 函 数 ， 同 样 ， 可 由 基本 全 点 形 函 数 经 过 收缩 
和 移动 来 得 到 ， 但 是 要 删除 超出 三 角形 单元 域外 L2 部 分 的 形 函 数 ， 所 以 ， 边 界 节 点 是 二 分 
节点 ;对 于 三 角形 单元 三 顶 节 点 〈T) 的 形 函数 ， 也 可 由 基本 全 点 形 函 数 经 过 收缩 和 移动 
来 得 到 ,但 是 要 删除 超出 三 角形 单元 域外 的 六 边 形 中 其 它 五 部 分 的 形 函 数 ， 所 以 , — DU SA 
都 是 分 节点 。 


图 4， 当 分 辩 率 RL=2X5 (n=3) 时 ， 三 角形 单元 的 节点 分 布 
Fig.4.The node partition for a triangular element when the RL=2 X 5 (n=3) 


HEAT a, OCA REISE. FE LETC, > Pras > 
Penns > Duas EAER n AAT PBS BOE BU "FJ ea A EAPC, PTCA SEPA DR 
I RE RT Fe NEETER, ENV, |) Ay f e AAP RS 
分 辨 率 分 析 Cmulti-resolution analysis, MRA), BẸ 


V,:= span(W, :ieZ) (12) 
E ur 


1 


其 中 ，2 WERSI V Os RL (11) (12) 的 位 移 子 空间 ， 


je, ERRA R= (n) (n2) 的 位 移 也 空间 中 ， 板 过 单元 内 任意 一 点 的 位 


移 可 写 为 ; 
u, 
v = Pra’ (13) 
w, 

式 中 

a | Diis 2 as a [uv E Ors m" [thane Mans Ws Oo Oo || 


Ue Vics Wey > orgs wawa (rst) En) ais y. z 方向 的 位 移 和 绕 x. 
n 


I n 
y 轴 的 转角 。 

由 此 可 以 看 出 ， 单 元 内 节点 是 规则 均匀 分 布 的 ， 其 数量 和 位 置 可 由 分 辩 率 常数 来 自 | 
调节 ; 一旦 确定 了 单元 分 辩 率 常数 也 就 规定 了 节点 总 数 及 其 具体 位 置 。 当 式 〈13 ) 中 的 分 辨 
率 常数 RL=1x3 (n=1) 时 ， 该 式 就 退化 为 传统 三 角形 板 壳 单元 的 位 移 场 ， 如 式 (4) Ara, 
因此 ， 传 统 三 角形 板 壳 单 元 是 单 分 辨 率 单元 ， 为 本 文 多 分 辩 率 三 角形 板 壳 单元 的 一 个 特例 。 


3 多 分 辩 率 三 角形 板 壳 单 元 


{EAP HON R= (n+ 1)(n + 2) 的 位 移 子 空间 中 ， 角形 板 壳 单 元 总 势能 可 写 为 ; 


nV,)=5 [fll Palle], ddv- ff, (p+ p.i) tna - Y [Ps Qno) Dag (4591) 


] 
+z ff, [xh [Ple], ddv- || qwidxdv- > Ow i) 
Pos p, DHA x, 了 轴 的 分 布 力 ; RAP P, 分 别 为 单元 上 某 处 的 沿 x, 了 轴 的 集中 
力 ; Q, 为 单元 上 某 处 的 横向 集中 荷载 ; q 为 横向 分 布 力 。 还 有 


Ou, ew 
x ex? 
ae} =] 98 |b, [e] =| Č |, 
n Oy n oy? 
Ou, | Ov, 5 we 
Oy Ox L OxOy 
l gu 0 [1 u 0 
[L]-G|» 1 9 [Blaa 1 (6 
0 0 (1-4)/2 |o 0 (1-4)/2 s 
EF 


C, = EC , E isthe FAMERS; 砷 板 壳 单 元 的 厚度 ， 4 为 泊 松 比 。 
-u 


K = [BM B7, B7 a: 


(15) 
[«], | 5 00* Bt BJ 
ð ap .| 
Pl 0 = S 
EX 中 B. = x 09 PA > a, = lu 00 sV 00 7" us Vis UC Un Vn ] 
0 n,rs n,rs 
Oy Ox 


e 


a, E [ woo d Oo 5 Ayo un Ws ^ 0, 5 0 "W, 0 0 li 


yrs ^ nn? "xnn? ^" ynn 


B? = 9$, C bsa uL i 
rs Ox? 5 y? E OxOy , 


将 式 (15)、(13) 代入 式 (14)， 并 经 整理 后 可 得 到 : 


II, (V,)= ua laK a — at? f* -a F° 


5 an Kia, 一 ax f, (16) 
Jub, K; A TERURERE RE: fy, REE p ASARI Fr AGRO netta 
——— 即 AT, (V,) = 0， 可 得 板 壳 单元 的 平衡 方程 : 

Keae =fe +F° (17) 
式 中 


k” : k” ; k” 
L 00 3 


rS nn 
aN 中 
rs 
k. = > k Js 
č P 
d-sļ|<1 
rs A = 
k. =k a, =9, when 


其 中 ，k ,为 节点 Ce d) 与 节点 (x,s) 的 相关 刚度 矩阵 ， 其 可 写 为 


= [[8., EB’ dxdy (18) 


D 


k. rs 


tz. = ([[*,.. ] adxay 
(19) 
Fe => [6,, (n -r.ny,-s)] P 


i 


式 中 N ZJ BEM; P 为 集中 荷载 癌 量 。 


c 4 ”转换 矩阵 


通过 上 述 步 又 , 可 以 得 到 多 分 辨 率 三 角形 板 壳 单元 的 刚度 与 等 效 节点 力 , 为 进行 结构 计 
FL, 需要 将 局 部 坐标 系 下 单元 节点 方向 ,转换 到 整体 坐标 系 下 单元 节点 方向 , 单 刚 与 等 效 节 
点 力 转换 到 整体 坐标 系 下 ， 此 转换 可 以 通过 转换 矩阵 工 来 实现 ， 即 


K' =T’K‘T’ (20) 
fi =T f° (21) 
Fi = TZ F° (22) 
其 中 
| 入， 0 
m cos 6, 0 0 
T; = À; À. =| 0 cosÜ,. cos@y 
0 cosy cosé 
0 入 ， 


其 中 0 为 单元 局 部 坐标 轴 与 整体 坐标 轴 之 间 的 夹 角 。 
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5 结构 总 刚 和 节点 力 


转换 到 整体 坐标 系 下 , 得 到 整体 坐标 系 下 的 单元 边界 分 节点 , 各 个 单元 相应 共同 边界 上 
的 分 节点 相互 全 加 (定义 域 扩 展 )， 就 可 易 得 到 结构 的 总 刚 和 节点 力 : 


K,=>K; (23) 
F,-YE (24) 
f,=)> f (25) 


6 算 题 验证 


题 1[7]. 一 两 对 边 简 支 的 平行 四 边 形 斜 板 , 各 边 长 为 工 ， 斜 交角 为 60， 如 图 5 Bras, Rr 
FAE = 0.3， 其 上 作用 有 均 布 荷载 49， 求 该 板 中 心 处 的 挠 度 。 


图 5. RHR 
Fig. 5. Askew plate 

就 斜 板 中 心 的 搁 度 求解 问题 ,我 们 分 别 采用 本 文 的 多 分 辨 紊 和 传统 的 单 分 辩 率 三 角形 板 壳 单 
元 有 限 元 方法 及 区 间 三 次 B 样 条 小 波 有 限 元 法 CBSWD 来 求解 。 由 于 采用 了 全 节点 形 函 数 ， 
多 分 辨 率 三 角形 板 壳 单 元 计算 模型 为 整体 结构 模型 ， 即 只 需 两 个 单元 GDG@) 来 模拟 整个 斜 板 ， 
每 个 单元 的 分 辨 率 常数 分 别 取 5x9、7x13 和 9x17， 其 分 辩 率 分 别 为 L/8、L/12 和 L/16, 
相应 的 整体 斜 板 的 分 状 率 常数 分 别 为 9x9、13x13 和 17x17, 斜 板 的 分 辩 率 与 单元 的 相同 ， 
分 别 为 L8、L/12 和 L116， 如 图 6a， 根 据 分 辨 率 的 定义 ， 要 求 单元 节点 均匀 分 布 ; 整体 斜 
板 节点 等 距 布置 , 因此 ， 一 旦 确定 了 分 辨 率 常 数 ， 那 么 ， 节 点 数量 和 位 置 也 就 确定 了 ; 然而 ， 
由 于 传统 有 限 元 采用 了 分 节点 形 函 数 ， 传 统 的 三 角形 板 壳 单 元 建立 的 计算 模型 是 离散 模型 ， 
其 离散 三 角形 网 格 划 分 取 对 应 的 8x8、12x12 和 16x16， 如 图 7b， 显 然 ， 这 里 的 网 格 节点 
分 布 形式 没有 均匀 性 的 要 求 ， 其 位 置 是 可 以 变动 的 ; 而 区 间 三 次 了 B 样 条 小 波 有 限 元 法 采用 
了 样 条 尺度 函数 矢量 积 的 方法 构建 了 复杂 的 全 节点 形 函 数 ， 所 以 该 模型 也 是 整体 结构 模型 。 
其 计算 结果 列 于 表 1.。 由 表 1. 看 出 ,我 们 可 以 通过 调节 多 分 辩 率 单元 分 辩 率 的 大 小 来 提高 计 
算 的 清晰 度 ， 此 单元 相当 于 一 张 图 片 ， 节 点 等 同 于 图 片 的 像素 ， 而 传统 单 分 辨 率 单 元 是 通过 
增加 网 格 的 密度 来 提高 计算 精度 , 而 精度 与 清晰 度 相当 , 但 是 与 传统 单 分 状 率 单元 有 限 元 法 
相 比 较 ， 多 分 辨 率 单 元 法 有 如 下 优势 : 多 分 辨 率 单元 模型 中 的 分 辨 紊 具有 明确 的 数学 意义 ， 
其 取 值 具 有 坚实 的 数学 基 ， 而 传统 模型 的 网 格 划 分 主要 依靠 经 验 ,， 无 数学 基 ， 无 任何 数学 意 
义 。2) 多 分 辨 率 整体 模型 中 只 有 两 个 单元 ， 因 此 ， 计 算 中 矢量 转换 次 数 少 〈 一 节点 转换 一 
次 )， 而 传统 离散 模型 中 单元 数量 多 ， 矢 量 转换 次 数 众 多 〈 一 节点 转换 六 次 )， 计 算 效 率 低 。 
比较 区 间 三 次 B 样 条 小 波 有 限 元 法 分 析 结 果 的 精度 ， 多 分 辨 率 单 元 法 的 计算 清晰 度 与 其 相 
当 ， 但 多 分 辩 率 单元 法 有 如 下 特长 D 


F 
— 
RRR 


(a) 多 分 辩 率 整体 模型 (b) 单 分 辨 率 有 限 元 离散 模型 
(a) The multiresolution integrated model (b) The monoresolution discretized model 
图 6. 斜 板 的 计算 模型 
Fig. 6 The computational model for the skew plate 


多 分 辨 率 单元 模型 中 的 分 辨 率 具 有 明确 的 数学 意义 , 其 取 值 具有 坚实 的 数学 基 , 而 区 间 三 次 
B 样 条 函数 节点 数量 的 选取 主要 依靠 经 验 ， 也 无 数学 基 ， 无 任何 数学 意义 。2) 全 节点 形 函 
数 表达 式 简洁 且 有 具有 Kronecker delta 特性 ,便于 处 理 边界 条 件 ， 从 而 提高 了 计算 效率 ， 而 矢 
量 积 构造 的 全 节点 形 函 数 表达 式 较为 复杂 ， 且 有 人 负 值 ， 不 利于 快速 数值 积分 ， 再 者 ， 它 没有 
Kronecker delta 特性 ， 因 此 需要 对 边界 条 件 进行 特殊 处 理 ， 从 而 大 大 地 影响 了 计算 速度 。 


表 1. 斜 板 中 心 的 挠 度 (w/ qL* /100Do) 
Table.1. the center point deflection (w/ qL^ /100Do) 


单元 类 型 
本 文 单元 传统 单元 TER 
(RL/ 单 个 单元 ，RL / 斜 板 ) (网 格 ) 
5x9, 9X9 8x8 0.7920 
7X13, 13x13 12x 12 0.7937 
9x17, 17X17 16X 16 0.7930 
BSWI ( He, 2006)[7] 0.7925 
Analytical (Hota, 1988)[7] 0.7945 


题 2[8]. 一 半径 为 x 厚度 为 ! 周 边 简 支 或 固 文 的 圆 形 薄 板 ， 承 受 数值 为 q 的 横向 均 布 压力 ， 
其 材料 弹性 模 量 为 妃 ， 泊 松 比 为 4&=0.3， 求 解 该 板 圆心 处 的 位 移 和 弯 矩 。 
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a A.multiresolution integrated model 


针对 求解 板 圆 心 处 位 移 和 弯 和 窍 的 问题 ， 利 有 


(a) 多 分 辩 率 整体 模型 


(b) “ 单 分 辨 率 有 限 元 离散 模型 


图 7. 14 圆 板 的 计算 模型 


Fig. 7. The computational model for the 1/4 circular slab 


进行 结构 计算 。 就 具体 的 计算 模型 ， 我 们 可 分 别 采 月 
统 的 单 分 辨 率 有 限 元 离散 模型 ， 如 图 7. 5)。 多 分 辨 率 整 体 模型 有 两 个 三 角形 单元 G@， 其 


b Amonoresolution discretized model 


结构 和 荷载 的 对 称 性 ,可 采用 四 分 之 一 圆 板 
多 分 辨 紊 整体 模型 ， 如 图 7. Ca), Tf 


每 个 单元 的 分 辨 率 常数 分 别 取 2x5、4x7; 传统 有 限 元 离散 模型 网 格 相 应 地 取 3x6、7x13。 


其 固 支 (SS) mfi (SC) 的 计算 结果 列 于 表 2。1 
周 节 分 辨 率 来 提高 计算 的 清晰 度 , 均 布 的 节点 数量 越 多 , 清晰 度 越 高 ;两 单元 
f 接 成 全 节点 ,从 而 完成 整体 结构 的 计算 。 与 传统 有 限 元 离散 法 相 比 


单元 法 可 通过 i 
可 在 交界 分 节点 处 进行 搓 


表 2 可 以 看 出 ， 多 分 辨 率 三 角形 板 壳 


较 ， 多 分 辨 率 三 角形 板 吝 单元 法 使 用 单元 数量 少 很 多 ， 因此， 矢量 转换 次 数 大 幅度 减少 ， 从 
而 大 大 地 提高 了 计算 速度 ; 另外 , 由 于 有 数学 基 的 支撑 , 用 分 辨 率 来 确定 节点 的 数量 和 位 置 ， 


以 及 描述 计算 的 精确 性 , 与 传统 的 网 格 划分 和 精度 描述 相 比 较 , 多 分 辨 率 三 角形 板 过 单元 法 


显得 更 加 科学 和 高 效 。 


R2 ARKE ORENSE 


Table.2.the Central deflection and Bending moment for the circular slab 


Central deflection (a x q/*/100 Do) Bending moment (8 xq/?/10) 
RL/elem | Mesh SS SC SC (central) 
Multi | Mono Multi | Mono Multi | Mono 
2x5 3X6 0.0145 | 0.0145 0.0638 0.0638 0.2103 0.2103 
4x7 7x13 | 0.0153 | 0.0153 0.0637 0.0637 0.2073 0.2073 
Analytical[8] 0.0156 0.0637 0.20625 


题 3[9]. 如 图 8 所 示 折 板 ， 弹 性 模 量 E-2.06x 10°Mpa, WAHE u = 0.3 ， 折 板 厚 度 1m, 


折 板 为 方形 板 ， 边 长 上 =50m。 受 均匀 外 压 g=1KN/m2 TER 


各 边 简 文 。 求 解 沿 图 8 中 AB 和 CD 边 的 位 移 。 


图 8， 折 板结 构 


Fig 8. A folded plate structure under the uniform loading q 


月 ， 边 界 条 件 为 AB、EF 自由 ， 其 余 
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a 多 分 辨 率 整 体 计算 模 型 


a Multiresolution integrated computational model 


本 题 分 别 采 用 本 文 的 多 分 辩 率 三 角 
单元 求解 。 为 了 满足 单元 内 布点 ] 
分 别 取 中 点 C、D， 由 C、D ERE ZH AO, XE, 
正方 形 区 域 , 每 个 正方 形 区 域 再 以 三 角形 划分 ,这样 ， 折 板 就 可 采 月 


b 单 分 辩 率 高 散 计 算 模型 


b Monoresolution discretized computational 


图 9 折 板 结构 计算 模型 


Fig 9. The folded plate structural computational model 


的 匀 ， 整 体 结构 布点 等 


单元 和 传统 有 限 元 中 的 单 分 辩 率 三 角形 板 壳 

E 的 分 辨 率 要 求 ， 在 该 折 板 两 底 边 上 
就 将 整个 折 板 分 成 四 个 同等 大 小 的 
H 8 个 多 分 辨 率 三 角形 板 


过 单元 DOS@@O@OB 组 成 的 多 分 辩 率 整体 计算 模型 ， 如 图 9a 所 示 ， 每 个 单元 的 分 辨 率 
常数 为 3x7 (n=5)， 而 整个 折 板 结构 的 分 辨 率 常数 是 11x 11， 整 体 结构 节点 均匀 且 等 距 分 


布 ， 因 此 ， 其 分 辨 率 为 Sm。 由 于 该 模型 域 中 节点 是 全 节点 ， 所 以 ， 它 属于 无 网 格 整体 计 入 
模型 ， 其 计算 结果 对 应 图 10a、 图 10b 中 的 case A HHA; ë 
传统 单 分 辩 率 三 角形 板 沉 单元 构成 ， 如 
它 属于 网 格 离散 计算 模型 ， 其 计算 结果 对 应 图 10a、 图 10b 中 的 case B 曲线 ; 图 10 中 的 第 
三 条 曲线 是 文献 [9] 采 用 区 间 三 次 B 样 条 小 波 有 限 元 分 析 的 结果 ， 其 计算 模型 是 | 
建 的 全 节点 形 函 数组 合成 的 无 网 格 整体 模型 ， 由 了 
性 ， 所 以 该 数学 模型 在 处 理 边界 条 伯 


H 
m 


图 9b 所 示 ， 由 于 该 模型 域 中 节点 


该 模型 的 计算 效率 ,模型 图 可 以 参见 文献 [9]。 晶 
元 法 既 有 有 限 元 法 的 特长 又 有 无 网 格 法 的 优点 , 它 既 可 以 根据 结构 几何 图 形 和 受 力 特性 
分 区 域 单元 ,又 可 以 通过 调节 其 单元 分 辩 率 的 大 小 来 调整 
计算 清晰 度 。 虽 然 多 分 辨 率 单元 模型 与 传统 的 和 


单元 节点 的 数量 ,从 而 调节 
元 模型 比较 ， 其 计算 清晰 度 相 当 ， 


F- 


是 由 10x 10 个 
节点 ， 所 以 ， 


REFR 
六 此 全 市 点 形 函 数 没 有 Kronecker delta 特 

上 方面 ， 需 要 转换 和 矩阵， 显得 比较 复杂 ， 因 此 大 大 降低 了 
HE 104,b 可 以 看 出 ， 三 角形 多 分 辨 率 板 壳 单 
ERRI 
元 的 


但 多 分 辩 率 模型 比较 传统 的 单 分 辨 率 模型 具有 以 下 优势 : CD 多 分 辨 率 模型 中 的 分 辨 率 具 有 


明确 的 数学 意义 ， 其 取 值 具有 坚实 的 数学 基 ， 而 传统 模型 的 网 格 划 分 3 


要 依靠 经 验 ， 无 数学 


基 , 无 任何 数学 意义 。(2) 由 于 域内 是 全 节点 ,多 分 辨 率 整 体 模型 中 单元 节点 刚度 与 等 效力 


可 以 直接 自动 生成 而 传统 离散 模型 的 单元 节点 为 分 节点 , 所 以 , 其 每 个 节点 刚度 与 等 效力 都 


需 进行 人 工 车 加 组 合 ; (3 ) 与 传统 单元 的 碎片 化 模型 相 比 较 , 多 分 辩 率 单元 的 整体 化 模型 中 ， 


单元 数量 少 ， 所 以 矢量 转换 次 数 少 ， 
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a. ”沿边 界线 AB 的 挠 曲 曲线 b，” 沿 中 线 线 CD 的 挠 曲 曲线 
a. The deflection curve along boundary line AB b. The deflection curve along midline CD 
图 10， 折 板结 构 位 移 反应 曲线 
Fig. 10. The folded plate structural displacement response curve 
(o 7 结论 
本 文 提 出 了 多 分 辨 率 三 角形 板 壳 单 元 法 ， 和 传统 板 壳 三 角形 单元 有 限 元 法 相 比较 ,本文 
e 方法 有 如 下 特点 : 
- D 提出 了 三 角形 板 壳 单 元 的 分 节点 和 全 节点 的 概念 , 并 基于 传统 三 角形 板 壳 单 元 的 分 
© 节点 形 函 数 ， 提 出 一 种 函数 三 角形 定义 域 的 延 拓 方 法 ， 构 建 了 基本 全 节点 形 函 数 ， 
=s 该 函数 具有 Kronecker delta 特性 , 便于 边界 条 件 的 处 理 , 同时 基本 全 节点 形 函 数 的 
= 连续 性 为 传统 有 限 元 法 中 节点 刚度 可 由 共同 节点 的 单元 刚度 县 加 得 到 的 处 理 方法 
提供 了 坚实 的 理论 依据 。 


= 2) 采用 基本 全 节点 函数 的 收缩 和 平移 , 构造 了 基 函 数 , 也 就 是 多 分 辨 分 析 的 简洁 而 坚 
实 的 数学 基 , 并 由 此 形成 了 多 分 辩 率 单元 的 位 移 子 空间 序列 , 并 引入 了 分 辨 率 常 
RL， 从 而 建立 了 理性 的 多 分 辨 分 析 概念 ， 实 现 了 结构 理性 化 的 多 分 辨 数值 分 析 。 

3) 传统 三 角形 板 壳 单 元 是 单 分 辩 率 单元 ， 为 本 文 多 分 辩 率 三 角形 板 壳 单 元 的 一 个 特 


4) 分 辩 率 定义 了 区 域内 节点 的 等 距 均 匀 分 布 。 结 构 数 值 分 析 清 晰 度 的 高 低 是 由 分 辩 率 
大 小 决定 的 而 非 网 格 的 密 丽 ;高 清晰 度 单元 的 节点 数 多 于 低 清晰 度 单 元 ; 

5) 本 文 三 角形 多 分 辩 率 板 壳 单 元 方法 计算 效率 高 于 相应 的 其 它 有 网 格 或 无 网 格 板 壳 
单元 有 限 元 法 ; 多 分 辩 率 单元 法 以 清晰 度 取 代 了 其 它 方法 的 精度 ， 以 分 辩 率 代替 了 
网 格 或 条 带 , 以 理性 多 分 辩 分 析 蔡 换 了 非 理性 多 分 辩 分 析 , 以 全 节点 替代 了 分 节点 ， 
从 此 ， 结 构 的 计算 模型 实现 了 从 碎片 化 《离散 化 ) 到 整体 化 的 转变 ; 

6)” 随 着 多 分 辨 率 单元 法 4001.2.3] 的 提出 ， 理 性 多 分 辨 分 析 将 会 在 工程 实际 中 得 到 广泛 
的 应 用 。 


N 


SS 


参考 文献 : 

[1] ZIENKIEWICZ.O.C.,TAYLOR.R.L. The Finite Element Method[M]. seventh ed., Butterworth 
-Heihemann,London. 2013 

[2] BELYTSCHKO T., LU Y. Y., GU L. Element - free Galerkin methods[J], International Journal for 
Numerical Methods in Engineering. 1994, 37: 229-256. 

[3] KRYSL.P, BELYTSCHKO.T, SFLIB: A library to compute the element free Galerkin shape 
functions[ J], Computer Methods in Applied Mechanics and Engineering, 2001, 190(15-17): 2181-2205 

[4] KO.J, KURDILA A J, PILANT M S. A class of finite element methods based on orthonormal, compactly 
supported wavelet [J]. Computational Mechanics.1995,16:235—244. 

[5] HUGHES TJ.R, COTTRELL J.A., BAZILEVS Y , Isogeometric analysis: CAD, finite elements, 
NURBS, exact geometry and mesh refinement[]]. Computer Methods in Applied Mechanics and. 
Engineering. 2005,194(39-41): 4135-4195. 

[6] TOSHNIWAL D, SPELEERS H, HIEMSTRA R. R, HUGHES T. J. R, Multi-degree smooth polar splines: 
A framework for geometric modeling and isogeometric analysis[J], Computer Methods in Applied 

Mechanics and Engineering, 2017,316(Apr), 1005-1061. 

[7] HOTA V.S., RAO GANGA., CHAUDHARY V.K., Analysis of skew and triangular plates in 
bending[ J]. Comp & Struct, 1988.28, 223-235. 

[8] ERE, Fh 完全 无 闭锁 厚薄 板 通 用 弯 角形 单元 [中 .南通 大 学 学 报 ， 
2003,2(2):5-11.(WANG-.Haixia, CEN-Song. Locking-free thick plate bending triangular 
element[J]. Journal of Nantong University, 2003, 2(2) :2-5. (in Chinese)) 


[9] 何 正 嘉 , 陈 雪 峰 , 李兵 , 等 . 小 波 有 限 元 理论 及 其 应 用 [M]. 北京 : 科学 出 版 社 , 2006.(HE-Zhengjia， 
CHEN-Xuefeng, LI- Bing, etc. Theory and engineering application of wavelet finite element 
method[M]. Beijing:Science Press 2006 (in Chinese)) 

[10] XIA Y M, LIU Y X, CHEN S L, et al. A rectangular shell element formulation with a new 
multi-resolution analysis[J]. Acta Mechnica Solida Sinica,2014 27(6): 612-625. 

[11] XIA Y M, CHEN S L. A hexahedron element formulation with a new multiresolution analysis[J]. 
Science China Physics, Mechanics, Astronomy, 2015,58(1): 014601-10. 

[12] 夏 逸 鸣 ， 唐 敢 ， 江 世 永 . 多 分 辨 率 Timoshenko 梁 单 元 []. 固体 力学 学 报 ，2014, 35 (1): 
57-62.(XIA-Yiming, TANG-Gan, JIANG Shiyong. Multiresolution Timoshenko beam element[J]. 
Chinese Journal of Solid Mechanics.2014,35(1):57-62.(in Chinese) 

[3] RA. £ 4 HER DO se Sñ ú 8 BR IE]. 中 国 科学 〈 物 理学 、 力 学 、 天 文学 ) , 2017, 
47(8 ) :084601-12 (XIA-Yiming. Multiresolution finite element based on a new quadrilateral shell 
element[J]. Science China-Phy Mech Astron,2017,47(8):084601-12.(in Chinese)) 


17 


A multiresolution triangular shell element method 


XIA-Yiming 
(1. School of civil aviation, Nanjing University of Aeronautics and Astronautics, Nanjing 211106, P.R.China) 


Abstract: Initially in accordance with domain scope of a traditional node shape function, a novel concept of a 
split/full node is presented. Based on a split node shape function for a triangular shell element, the basic full 
node shape function is constructed by extending the split node shape function to other three quadrants around 
the coordinate zero point and a displacement subspace sequence is constituted out of scaling and shifting of the 
basic full node shape function on the element domain, which brings about a simple, clear and rigorous 
mathematical basis for the multiresolution analysis (MRA).After that, the multiresolution triangular shell element 
and method is formulated through applying the minimum potential energy principle. Hence, the MRA concept is 
established and the resolution level (RL) for the element is introduced. Finally, through a numerical example, the 
following conclusions is drawn: The RL defines an equally-spaced uniform node distribution layout. The traditional 
3-node triangular shell element and method is essentially a mono-resolution one and also a special case of the 
proposed element and method. Due to the full nodes in the element domain, the structural model for a 
numerical analysis by the proposed multiresolution element is an integrated one by RL, not a discretized one by 
mesh just as by the conventional monoresolution element. The proposed finite element method is a rational MRA 
method and overperforms all those irrational MRA methods in numerical analysis efficiency, also can consolidate 
those irrational MRA methods. The clarity of a structural analysis is defined by the RL, not mesh. The continuity of 
the full node shape functions unveils the secretes behind artificial assemblage of the corresponding items of the 
stiffness，mass matrices and the equivalent node loading vectors at common nodes by the conventional 


quadrilateral shell element. 


Key words: triangular shell element; split node; full node; displacement subspace sequence; rational 


multiresolution analysis; resolution level (RL) 
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